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In the last decade, transition metal doped ZnO has been intensively investigated as a route to 
room temperature diluted magnetic semiconductors (DMS). However the origin for the reported 
ferromagnetism in ZnO based DMS remains questionable. Possible options are diluted magnetic 
semiconductors, spinodal decomposition or secondary phases. In order to clarify this question, we 
have performed a thorough characterization of the structural and magnetic properties of Co and 
Ni implanted ZnO single crystals. Our measurements reveal that Co or Ni nanocrystals (NCs) are 
the major contribution of the measured ferromagnetism. Already in the as-implanted samples, Co 
or Ni NCs have formed, and they exhibit superparamagnetic properties. The Co or Ni NCs are 
crystallographically oriented with respect to the ZnO matrix. Their magnetic properties, e.g. the 
anisotropy and the superparamagnetic blocking temperature can be tuned by annealing. We discuss 
the magnetic anisotropy of Ni NCs embedded in ZnO concerning the strain anisotropy. 



I. INTRODUCTION 

Recently diluted magnetic semiconductors (DMS) are 
under intensive investigation due to their potential ap- 
plications in spintronics. In DMS materials, transition 
or rare-earth metal ions are substituted onto cation sites 
and are coupled with free carriers to yield ferromagnetism 
via indirect interaction. Mn-doped InAs [1 and GaAs [2] 
are the best studied DMS materials. Conventional III- 
V semiconductors are widely used for high-speed elec- 
tronic and optoelectronic devices. The discovery of 
hole- mediated ferromagnetism in (Ga,Mn)As opened the 
way to integrate spin-based phenomena into mainstream 
microelectronics and optoelectronics as well as taking 
advantage of already established fabrication processes. 
Electrically controlled spintronic devices based on GaM- 
nAs and InMnAs have been successfully designed and 
tested at low temperatures, e.g., a spin-LED using GaM- 
nAs as a spin injector [3]. However the highest Curie 
temperature (Tc) reported in (Ga,Mn)As is ^170 K [5, 
which is far below room temperature and limits its reg- 
ular application. 

In 2000, Dietl et al. [5 used a mean field theory to 
estimate the ordering temperature Tc of ferromagnetic 
semiconductors, and they predict that room-temperature 
ferromagnetism can be created by substituting Mn ions 
in p-type wide-band gap semiconductors such as GaN 
and ZnO. Sato et al. calculated the properties of n-type 
ZnO doped with M TM ions (V, Cr, Mn, Fe, Co, and Ni) 
[6j. The ferromagnetic state, with a Tc of around 2000 
K, is predicted to be favourable for V, Cr, Fe, Co, and Ni 
in ZnO while Mn-doped ZnO is predicted to be antiferro- 
magnetic. These predictions largely boosted intensive ex- 
perimental activities on transition metal doped GaN and 
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ZnO. A large number of research groups have reported 
the experimental observation of ferromagnetism in TM 
(from Sc to Ni) doped ZnO H El [1 [lOl [III [HI [13 [H [15] 
fabricated by various methods including ion implanta- 
tion. For a comprehensive review, see Ref. p!6| [T7| [TS]. 
However in these reports the magnetic properties us- 
ing the same dopant vary considerably. E.g. the sat- 
uration moment and Curie temperature for Mn doped 
ZnO ranges from O.OTS/i^/Mn, 400 K [9 to O.lT/i^/Mn, 
30-45 K [12], respectively. In contrast to these publica- 
tions, other groups reported the observation of antiferro- 
magnetism [19 , [20l [21], spin-glass behavior [22 , 23 , and 
paramagnetism [20l [21 [351 [26] in TM-doped ZnO. Re- 
cently it was also found that nanoscale precipitates can 
contribute to the ferromagnetic properties substantially 
[27l[2i[29l[30l|3ll|3a|33|34]. 

Publications claiming the intrinsic ferromagnetism in 
TM doped ZnO are often based solely on magnetization 
measurements using high sensitivity superconducting 
quantum interference device magnetometery (SQUID) 
and structure characterizations using lab-equipped x-ray 
diffraction. The latter has been demonstrated to be not 
sensitive enough to detect nanoscale precipitates [28^. 
Nowadays the research community puts more effort to 
judge if TM ions in semiconductors are homogenously 
distributed, and if the ferromagnetism is intrinsic [35 , 36]. 
Anomalous Hall effect (AHE), that verifies charge car- 
rier participation in the magnetic order, has been sug- 
gested to judge about carrier-mediated ferromagnetism 
arising from DMS. However a recent study presents the 
co-occurrence of superparamagnetism and AHE in Co 
doped Ti02 films [37]. Another criterium is the mag- 
netic anisotropy [TH [38l |39] as a signature of the in- 
trinsic ferromagnetism. The controversy in the magnetic 
properties of ZnO-based DMS, as stated above, might 
partially be due to the insufficient characterization of 
the samples [^ • Particularly, a careful correla- 
tion between structure and magnetism should be estab- 
lished by sophisticated methods. Synchrotron radiation 
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based x-ray diflFraction (SR-XRD) is a powerful tool to 
detect small precipitates, e.g. metallic TM nanocrystals 
in ZnO [28]. On the other hand, element selective mea- 
surements of the magnetic properties, e.g. XMCD [43J, 
and Mossbauer spectroscopy [28, 32^, address the origin 
of the measured magnetism directly. 

The present work is devoted to a comprehensive in- 
vestigation of Co and Ni doped ZnO. By SR-XRD and 
SQUID magnetometry, we correlate the structural and 
magnetic properties. We attempt to answer the follow- 
ing questions concerning the presence of ferromagnetic 
precipitates. 

(1) How do the crystalline precipitates orient with re- 
spect to the host matrix? 

(2) Which techniques are suitable to detect these pre- 
cipitates? 

(3) Are these precipitates the major contribution of 
the measured ferromagnetism, or is there another source 
which contributes? 

(5) Can these nanoscale precipitates exhibit magnetic 
properties, e.^. magnetic anisotropy, concerning the crite- 
ria above mentioned? 

(6) How do these precipitates behave upon thermal 
annealing? 

The paper is organized as follows. In section |TT] all 
the experimental methods employed will be described. 
Then the results will be presented in two sections. In 
section HI we will focus on the as-implanted samples. 



and discuss the orientation, the superparamagnetism, the 
magnetic anisotropy of Co and Ni NCs. In section IV we 



will describe the structure and magnetism evolution due 
to high vacuum annealing. Finally in section|V|we discuss 
the origin of the magnetic- anisotropy for the oriented Co, 
Ni NCs system, and the possible formation of Co/CoO 
and Ni/NiO core/shell structures upon annealing at 923 



K. This paper is concluded in section VI 



II. EXPERIMENTS 



backscattering yield at channeling condition to that for 
a random beam incidence [45 . Therefore, Xmin 

labels 

the degree of lattice disorder upon implantation, and an 
amorphous sample shows a Xmin of 100%, while a perfect 
single crystal corresponds to a Xmin of 1-2%. 

Magnetic properties were measured with a supercon- 
ducting quantum interference device (SQUID, Quantum 
Design MPMS) magnetometery. The samples were mea- 
sured with the field along both the in- and out-of-plane 
direction. We studied both the temperature dependence 
of the magnetization at a constant field and the field de- 
pendence at a constant temperature. By SQUID, virgin 
ZnO is found to be purely diamagnetic with a susceptibil- 
ity of -2.65x10"^ emu/Oe-g. This background has been 
subtracted from the magnetic data. The temperature 
dependent magnetization measurement has been carried 
out in the following way. The sample was cooled in zero 
field from above room temperature to 5 K. Then a 50 Oe 
field was applied, and the zero field cooled magnetization 
curve (ZFC curve) was measured with increasing temper- 
ature from 5 to 350 K, after which the field cooled mag- 
netization curve (FC curve) was measured in the same 
field from 350 to 5 K with decreasing temperature. 

Structural analysis was performed both by synchrotron 
radiation x-ray diffraction (SR-XRD) and conventional 
XRD. SR-XRD was performed at the Rossendorf beam- 
line (BM20) at the ESRF with an x-ray wavelength of 
0.154 nm. Conventional XRD was performed with a 
Siemens D5005 equipped with a Cu-target source. In 
XRD measurement, we use 20-0 scans to identify crys- 
talline precipitates, and pole figures (azimuthal (/)-scan) 
for determining their crystallographical orientation. As 
a standard approach, for an XRD ^-scan one first tilts 
the sample by the angle of x from the sample surface 
(x is the angle between the diffraction plane of inter- 
est and the sample surface), and fixes the Bragg angle. 
Subsequently the spectrum is recorded during azimuthal 
rotation with respect to the sample normal. The pole 
figure is constructed by a series of 0-scans at different x- 



Commercial ZnO single crystals grown by the hy- 
drothermal method were implanted with Co or Ni ions 
at 623 K with a fiuence ranging from 0.8 x 10^^ cm~^ 
up to 8 X 10^^ cm~^. The implantation energy was 180 
keV, which resulted in a projected range of Rp = 89±29 
nm, and a maximum atomic concentration from ~1% to 
^10% (TRIM code [44 ). Thermal annealing was per- 
formed in a high vacuum (<10~^ mbar) furnace from 
823 K to 1073 K. 

The lattice damage induced by implantation was evalu- 
ated by Rutherford backscattering/channeling spectrom- 
etry (RBS/C). RBS/C spectra were collected with a col- 
limated 1.7 MeV IIe+ beam at a backscattering angle of 
170°. The sample was mounted on a three- axis goniome- 
ter with a precision of 0.01°. During channeling measure- 
ment, the sample was aligned to make the ZnO<0001> 
axis parallel to the impinging He+ beam. Xmin^ the 
channeling minimum yield in RBS/C, is the ratio of the 



III. AS-IMPLANTED SAMPLES 

A. ZnO lattice damage upon implantation 

RBS/C is used to check the lattice damage after im- 
plantation. Figure [l] shows RBS/C spectra for different 
fluences. The arrow labelled Zn indicates the energy for 
backscattering from surface Zn atoms. The implanted 
Co or Ni ions cannot be detected for the very low fiu- 
ence (0.8x10^^ cm~^, not shown). However, they are 
more pronounced as a hump in the random spectrum for 
a larger fiuence of 4x10^^ cm~^ and 8x10^^ cm~^ (not 
shown). The humps in the channeling spectra mainly 
originate from the lattice disordering due to implanta- 
tion. As expected, Xmin increases with increasing fiuence 
(see Table |l|. Note that the highest Co fiuence induced 
comparable lattice damage with the middle fiuence, and 
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FIG. 1: RBS random (ran.) and channeling (ch.) spectra, (a) 
Co implanted ZnO, and (b) Ni implanted ZnO (The fluence 
for Co and Ni ions is indicated). The yield of channeling 
spectra is progressively decreased with increasing fluence. 
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FIG. 2: SR-XRD 20-0 scan reveahng the existence of Co or Ni 
precipitates in (a) Co implanted ZnO, and (b) Ni implanted 
ZnO. 



less than the damage created by the same fluence of Ni. 
The reason could be a dopant specific self- annealing pro- 
cess [46^. RBS/C measurements also reveal that the ZnO 
is a hard material with respect to irradiation. The host 
material still partly remains in a crystalline state after ir- 
radiation with Co and Ni ions up to a fluence of 8x10^^ 
cm~^ {Xmin of 59% and 69%, respectively). 



B. Crystalographically oriented Co and Ni NCs 

SR-XRD is used to identify the precipitates in ZnO af- 
ter Co or Ni implantation. Figure |2] shows the XRD 20-0 
scans for all samples implanted with different fluences. 
At a low fluence (0.8x10^^ cm~^), no crystalline Co or Ni 
NCs could be detected. At a fluence of 4x10^^ cm~^ the 
hcp-Co(0002) (or fcc-Co(lll)) and Ni(lll) peak appear, 
respectively, and grow with increasing fluence. The full 
width at half maximum (FWHM) of the Co or Ni peak 
decreased with fluence, indicating the growth of the av- 
erage diameter of these NCs (table The crystallite size 
is calculated using the Scherrer formula ^47^, 



(i = 0.9A/(/3-cos(9), (1) 

where A is the wavelength of the x-ray, the bragg angle, 
and (3 the FWHM of 20 in radians. 

Note that there is only one peak for Co or Ni de- 
tectable, which indicates that the crystallites of Co 
or Ni are highly oriented with respect to the host 
matrix. The surface orientation is hcp-Co(0001)(or 
fcc-Co(lll))||ZnO(0001) and fcc-Ni(lll)||ZnO(0001), re- 
spectively. 

The Co hep structure only differs in the stacking from 
the fee one. Since the Bragg angles {0) for hcp-Co(0002) 
(l9=22.38°) and fcc-Co(lll) (l9=22.12°) are rather close 
to each other, it is difficult to assign these peaks in Figure 
|2|a) to hcp-Co or fcc-Co. A (/)-scan or pole figure on one 
of the diffraction planes not parallel with the sample sur- 
face (i.e. tilted by an angle x from sample surface) helps 
to identify hep or fcc-Co NCs and also reveals the crys- 
tallographical orientation relationship. By this approach, 
we find that only hcp-Co is present in the as-implanted 
samples. Figure [s] (a) and (b) show the pole figure of 
hcp-Co(lOll) and Ni(200), respectively. The radial co- 
ordinate is the angle (x) by which the surface is tilted 



FIG. 3: XRD Pole figure revealing the crystallographical ori- 
entation relationship between Co/Ni NCs and ZnO matrix, 
(a) Co(lOll) (in square) together with the tail of ZnO(1012) 
(in circle); The data points out of squares and circles are due 
to the background; (b) Ni(200) (in square) together with the 
tail of ZnO (101 2) (in circle). 
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FIG. 4: Hysteresis loops measured at 5 K for Co implanted 
ZnO with different fluences. Inset: Hysteresis loop measured 
at 300 K for the sample with the highest fluence. 



out of the diffraction plane. The azimuthal coordinate 
(^) is the angle of rotation about the surface normal. 
The pole figure shows poles of hcp-Co(lOll) at x~61.9°, 
and Ni(200) at x~54.8°, respectively. Both exhibit a 
sixfold symmetry. Since ZnO (1012) and hcp-Co(lOll) 
have similar Bragg angle, the poles of ZnO (1012) also 
show up at x~42.8° with much more intensities. The 
results are consistent with the theoretical Co(lOll) pole 
figure viewed along [0001], and Ni(200) pole figure viewed 
along [111] direction, respectively. Therefore, we can con- 
clude that these Co and Ni NCs are crystallographically 
oriented with respect to the ZnO matrix. The in-plane 
orientation relationship is hcp-Co[1010] ||ZnO[1010], and 
Ni[112] ||ZnO[1010], respectively. Due to the hexagonal 
structure of Co and sixfold symmetry of Ni viewed along 
[111] direction, it is not difficult to understand their crys- 
tallographical orientation onto hexagonal- ZnO. 

At this point, we have to remind the reader that in 
our former work on Ni implanted ZnO, the XRD mea- 
surement was performed on a conventional lab-equipped 
diffractometry (CXRD) [29]. CXRD reveals similar re- 
sults as SR-XRD (from the Ni fiuence of 4x10^^ cm-^ 
Ni NCs start to form). However, CXRD fails to detect 
Fe NCs in ZnO [28^, where Fe NCs are not well ori- 
ented like the case of Ni in ZnO. The peak intensity in 
XRD is proportional to the diffraction volume, i.e. the 
number of diffraction planes which are parallel to each 
other. For a 20-9 scan, the crystallographic orientation 
results in much more diffraction volume than the ran- 
dom orientation. Therefore, the high ordered orientation 
of NCs make them easier to be detected. As expected, 
Co NCs are also detectable in the CXRD measurement 
(not shown). 



C. Magnetic properties of Co and Ni NCs 

The magnetic properties of Co and Ni implanted ZnO 
were measured by SQUID magnetometery with the field 
parallel and perpendicular to the sample surface. Co im- 
planted samples exhibit a hard axis parallel to the sam- 
ple surface, while the easy axis is perpendicular to the 
surface. For Ni-implanted samples, the anisotropy direc- 
tions are vice versa. In this section, we investigate the 
superparamagnetism of the implanted ZnO samples and 
their magnetic anisotropy. 

1. Superparamagnetic Co and Ni NCs 

From the XRD results, we know that Co and Ni NCs 
have been formed in the as-implanted samples. For mag- 
netic nanoparticles, the formation of domain walls is en- 
ergetically unfavorable. Below a certain size (typically 
in the range of 15 to 30 nm depending on the mate- 
rial), the particle stays in a single-domain configuration. 
If the particle size is sufficiently small, above a partic- 
ular temperature (so-called blocking temperature, T^) 
thermal fiuctuations dominate and no preferred magne- 
tization direction can be defined. Such a system of su- 
perparamagnetic particles does not exhibit a hysteresis 
curve above T^; therefore the coercivity (He) and the 
remanence (M^^) are both zero. 

For a dc magnetization measurement in a small mag- 
netic field by SQUID, T^ is given by 

TB,S,md-^^^. (2) 

where K^ff is the anisotropy energy density, V the par- 
ticle volume, kB the Boltzmann constant [48]. For bulk 
crystals, Keff{V) is 5.0x10^ and 5.7x10^ Jm-^ for Co 
and Ni, respectively, at room temperature. 



5 



TABLE I: Structural and magnetic properties for Co and Ni-implanted ZnO with different fluence. Metallic Co/Ni fraction 
corresponds to the percentage of metallic Co/Ni compared with all implanted Co/Ni. 



Fluence 




Crystallite size 


Saturation magnetization^ 


Metallic 


Coercivity^ 


Is 


Orystallite size 


(cm ) 


(KBb/Cj 


(nmj (aKDJ 


[fJ^B/Co or /Ni) 


traction 


(Oe) 


(K) 


(nmj 


Co: 0.8x10^^ 


44% 














Co: 4x10^^ 


54% 


5 


0.29 (5 K) 


17% 


1400 (5 K) 


45 


4.3 


Co: 8x10^^ 


57% 


8 


0.44 (5 K) 


26% 


1400 (5 K) 


300 


8.1 


Ni: 0.8x10^^ 


45% 




0.05 (10 K) 


8% 


10 (10 K) 


< 5 


< 9.4 


Ni: 4x10^^ 


57% 


6 


0.16 (10 K) 


27% 


30 (10 K) 


16 


14 


Ni: 8x10^^ 


69% 


8 


0.22 (10 K) 


37% 


120 (10 K) 


70 


23 



^Refer to the easy axis at a saturation field of 10000 Oe for Co- 
ZnO, and 1500 Oe for Ni-ZnO, respectively. 
^Calculation from the average blocking temperature by Eq. I2I 



Phenomenologically there are two characteristic fea- 
tures in the temperature dependent magnetization of a 
nanoparticle system. One is the irreversibility of the 
magnetization in a small applied field {e.g. 50 Oe) after 
zero field cooling and field cooling (ZFC/FC) [48 . The 
other is the drastic drop of coercivity and remanence at 
a temperature close to or above [37]. 

Figure |4] shows the magnetization versus field rever- 
sal (M-H) measured at 5 K with the field applied per- 
pendicular to the sample surface (along ZnO [0001]). A 
hysteretic behavior is observed for the high-fluence im- 
planted samples. A saturation behavior is also observed 
at 300 K for the sample with the highest fluence (Figure 
[4] inset). However, neither coercivity nor remanence can 
be observed at 300 K. This is a strong indication for the 
superparamagnetism of a magnetic nanoparticle system. 
Knowing the formation of hcp-Co from XRD, it is reason- 
able to conclude that hcp-Co NCs are responsible for the 
magnetic behavior. For bulk hcp-Co crystals, the mag- 
netic moment is 1.7 /i^/Co at K. Assuming the same 
value for Co NCs, around 17% and 26% of implanted Co 
ions are in the metallic state for the fluence of 4x and 
8x10^^ cm~^ respectively. Similar results are observed 
for Ni implanted ZnO [29 . 

Note that the hysteresis loop for the fluence of 4x10^^ 
cm~^ exhibits two reversal steps. There is a very small 
kink when the field decreases from -2000 to -4000 Oe 
(and increases from 2000 to 4000 Oe). This wasp- waist 
shape of the loop is associated with magnetic phases with 
different coercivities jSTl [49] . 

Temperature dependent magnetization with H= 50 Oe 
was measured after ZFC and FC to confirm the super- 
paramagnetism (Figure [5|. For the magnetic samples, a 
distinct difference in ZFC/FC curves is observed. ZFC 
curves show a gradual increase (deblocking) at low tem- 
perature, and reach a broad peak, while FC curves con- 
tinue to increase with decreasing temperature. The broad 
peaks in ZFC curves are due to the size distribution of 
Co NCs. In this paper, the temperature at the maximum 
of the ZFC curve is taken as the average blocking tem- 



perature (later referred as T^). The ZFC/FC curves are 
general characteristics of magnetic nanoparticle systems 
with a broad size distribution [50[. T^ increases with 
the fluence, i.e. the size of nanoparticles. Table [T| lists 
the average size of Co NCs calculated by Eq. [l] (XRD) 
and Eq. [§] (SQUID). They are in a good agreement al- 
though there is also a size distribution in Co NCs [STj. 
The ZFC/FC magnetization was also measured for Ni 
implanted samples [29|. Comparing with Co, Ni has a 
much lower anisotropy energy density. For similar sizes 
of Ni NCs, the blocking temperature is therefore much 
lower than that of Co. Table |l] lists the average size of 
Ni NCs calculated by Eq. [l](XRD) and Eq. [2] (SQUID). 
Although the trend is the same for both calculations, the 
values from SQUID are much larger than that from XRD 
data. One reason could be the anisotropy energy density 
is underestimated by assuming the magneto-crystalline 
anisotropy constant. Another reason is that T^ is over- 
estimated by taking the temperature at the maximum of 
the ZFC curve. According to the calculation by Farle et 
al. [51 , and by Jacobsohn et al. [52 , the size dispersion 
for a given average particle diameter broadens the peak 
of ZFC curve and shifts it to higher temperatures. 



2. Magnetic anisotropy 

M-H loops were also measured for selective samples 
with Co fluence of 8x10^^ cm~^ with the field both per- 
pendicular and parallel to the sample surface. Figure 
|6ja) shows the comparison of the magnetization along 
ZnO[1010] and [0001] at 300 K. Figure |6|c) shows the 
orientation relationship between hcp-Co and ZnO, and 
the measurement geometry. Obviously Co [0001] is the 
easy axis, the same as a bulk hcp-Co crystal. The in- 
tersection of both curves gives an effective anisotropy 
field of 3000 Oe. At 5 K, the magnetic anisotropy is the 
same (not shown), and the coercivity of the easy axis is 
around 1400 Oe. The ratio between remanence and sat- 
uration magnetization is around 60%. This rather low 
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FIG. 5: Magnetization curves at 50 Oe after ZFC/FC for the 
Co implanted ZnO. With increasing fluence, the Co NCs is 
growing in size, resulting in a higher blocking temperature. 




FIG. 6: Hysteresis loops measured with the field along 
ZnO [0001] (out-of-plane) and [1010] (in-plane) for Co/Ni im- 
planted ZnO with the fluence of 8x10^^ cm~^ measured at 
300 K, (a) Co implanted ZnO; and (b) Ni implanted ZnO. 
(c) and (d) show the schematic geometry for magnetization 
measurements. 



remanence (below 100%) is due to the size distribution 
of nanomagnets. Very small nanomagnets behave su- 
perparamagnetically even at low temperatures, and only 
have field induced magnetization. This is a rather uni- 
versal feature for nanomagnets. For instance, epitaxial 
MnAs nanoclusters in GaAs, the remanence is also be- 
low 100% along the easy axis [53 . Figure [6]^b) shows the 
same measurement of Ni implanted ZnO, while (d) shows 
the orientation relationship between fcc-Ni and ZnO, and 
the measurement geometry. In contrast to bulk Ni where 
[111] is the easy axis, the easy axis is Ni[112] and the 
hard axis is Ni[lll]. Moreover, as shown in Figure [6]^b) , 
another in-plane direction Ni[110] is also an easy axis. 
Within the applied field, the magnetization curve along 
hard axis does not intersect with that along the easy axis. 
The effective anisotropy field is much larger than 1500 
Oe. That means that there are other contributions to 
the anisotropy dominating over the crystalline magnetic 
anisotropy. This will be discussed in section [V] 

In the work by Norton et al. [27 , epitaxial Co 
nanocrystals have been observed in Co implanted ZnO 
single crystals. The nanocrystal size is estimated to be 
~3.5 nm, which is below the superparamagnetic limit at 
room temperature. Therefore the ferromagnetism above 
300 K is very possible due to Co substitution onto the 
Zn site in the ZnO matrix. In our case, the measured 
superparamagnetism is well explained by the presence 
of Co and Ni nanocrystals. The formation of metallic 
nanocrystals already in the as-implanted state is due to 
the elevated implantation temperature (623 K), which 
facilities the precipitation process. Implantation at low 
temperatures {e.g. 253 K) prohibits precipitation, how- 
ever results in non-magnetism in our case, e.g. Fe im- 
planted ZnO [54 . 



IV. THE EFFECT OF POST- ANNEALING 

The post-annealing was performed in high vacuum 
with temperatures ranging from 823 K to 1073 K for 15 
min. The samples to be annealed were selected accord- 
ing to the SQUID measurement capability, namely the 
temperature range from 1.8 to 400 K. Due to the higher 
anisotropy energy of Co, T^ in ZFC curve of the high- 
est fluence is already 300 K in the as-implanted state 
(see Figure |5|. The annealing is expected to increase 
the size of Co NCs, and consequently increase T^ in the 
ZFC curve, which will exceed the temperature range of 
the SQUID magnetometery. Therefore, the sample im- 
planted with the middle Co fluence was chosen. For Ni 
the highest implantation fluence was chosen. 



A. Lattice recovery 

As shown in Figure [l] ion implantation substantially 
induces lattice damage of ZnO crystals. Here we show 
that this damage can be partially recovered by post- 
annealing. 

Figure [t] shows the RBS/C spectra for the sample im- 
planted with a Co fluence of 4 x 10^^ cm~^. The annealing 
temperatures are given in the figure. With increasing an- 
nealing temperature, the channeling spectra indicate that 
the lattice disorder of ZnO progressively decreases. After 
annealing to 1073 K, the channeling spectrum is almost 
comparable with the virgin ZnO. Similar RBS/C results 
are observed for Ni-implanted ZnO upon annealing (not 
shown) . 
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FIG. 7: RBS random (ran.) and channeling (ch.) spectra for 
Co implanted ZnO with a fluence of 4x10^^ cm~^ after ther- 
mal annealing at different temperatures. The lattice damage 
induced by implantation is progressively reduced by increas- 
ing annealing temperature. 



B. Evolution of structural properties 

Figure [sja) shows the development of Co NCs upon 
thermal annealing. The peak area and crystallite size 
calculated using the Scherrer formula [47 are compared 
in table T A broad scan (the inset of Figure [sja)) reveals 
only one peak from Co besides the ZnO peaks. An XRD 
0-scan has been used to distinguish between hep- and fcc- 
Co (Figure [9|. We find only hcp-Co in the as-implanted 
sample and the sample annealed at 923 K, while both fcc- 
and hcp-Co are present in the sample annealed at 823 K. 
The broad peak in Figure ga) (823 K ann.) is a su- 
perposition of hcp-Co(0002) and fcc-Co(lll). The crys- 
tallographical orientation relationship between Co NCs 
and ZnO is hcp-Co(OOOl) [1010] //ZnO(OOOl) [1010] //fcc- 
Co(lll)[112]. With this orientation, the 20-0 scans for 
hcp-Co(lOll) and fcc-Co(200) are expected with a skew 
geometry at one of the azimuthal position {e.g. at ^=0) 
as shown in Figure p^a) and (b), respectively. In skew 
geometry, the incident and the diffracted waves have the 
same angles to the surface, while the sample is titled with 
respect to its surface normal. By this configuration, a 
noncoplanar, its surface normal does not lie in the plane 
defined by the incident and the diffracted waves, can be 
measured [55 . Note that the peak area of Co in Fig- 
ure [Sj which is an approximate measure of the amount 
of Co NCs, increases drastically after 823 K annealing, 
while decreases after 923 K annealing. It is reasonable 
to attribute this change to the formation and disappear- 
ance of fcc-Co. The fcc-Co is probably oxidized to the 
amorphous CoO after 923 K annealing, while finally the 
majority of Co NCs are oxidized to amorphous state after 
annealing at 1073 K. 

For Ni implanted ZnO, the structure evolution upon 
annealing is similar as shown in Figure [sj^b) . The peak 
area and the crystallite size is listed in table III) The mild 
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FIG. 8: XRD 20-0 scans: (a) Co implanted ZnO crystals 
with different annealing temperature. The wide range XRD 
pattern for one of the samples (inset) reveals that no other 
crystalline phase (e.g. CoO) could be detected. The arrow 
points the peak shoulder coming from fcc-Co(lll) diffraction 
in the sample of 823 K ann. (b) Ni implanted ZnO crystals 
with different annealing temperature. The inset shows a com- 
parison between as-implanted and 823 K annealed samples. 



temperature annealing (823 K) only slightly increases the 
grain size of Ni. The annealing at 923 K drastically de- 
creases the peak area, while the grain size also decreases. 
Note that there is a significant shift in the Ni(lll) peak 
(labelled by dash lines), which indicates the relexation of 
lattice strain. After annealing at 1073 K, the majority of 
Ni NCs could be oxidized to amorphous state. 



C. Evolution of magnetic properties 

The structural phase transformation of Co NCs results 
in different magnetic properties as revealed by SQUID. 

^a) shows the ZFC/FC magnetization curves 
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Figure 

for all samples annealed at different temperatures. Ob- 
viously, except the sample annealed at 1073 K, the ZFC 
curves show a gradual increase (deblocking) at low tem- 
perature, and reach a maximum at a temperature of T^ 
(shown in table |llj), while FC curves continue to increase 
with decreasing temperature. No significant magnetiza- 



8 



TABLE II: Structural and magnetic properties for Co and Ni-implanted ZnO with different fluence. Metallic Co/Ni fraction 
corresponds to the percentage of metallic Co/Ni compared with all implanted Co/Ni. 



Sample 


Peak area 


Crystallite size 


Saturation magnetization^ 


Metalhc 


Coercivity^ 


Tb 


Crystallite size^ 




(AKDj 


(nmj (AKUj 


[fiB/^o or /Ni) 


traction 


(Oe) 


(K) 


(nm) 


Co: As-imp. 


380 


5 


0.29 (at 10000 Oe) 


17% 


1400 


45 


4.3 


Co: 823 K ann. 


8330 




0.36 (at 2000 Oe) 


21% 


250 


80 


5.2 


Co: 923 K ann. 


1400 


10 


0.32 (at 2000 Oe) 


19% 


450 


330 


8.3 


Co: 1073 K ann. 

















Ni: As-imp. 


7312 


8 


0.22 (at 2000 Oe) 


37% 


230 


70 


< 9.4 


Ni: 823 K ann. 


6590 


9 


0.22 (at 2000 Oe) 


37% 


200 


80 


14 


Ni: 923 K ann. 


750 


7 


0.18 (at 2000 Oe) 


30% 


220 






Ni: 1073 K ann. 


















^Refer to the easy axis magnetization at 5 K at a saturation field 
as indicated. 

^Calculation from the average blocking temperature by Eq. [3] 
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EIG. 9: XRD 0-scans for hcp-Co(lOll) (6'=23.78° and 
X'-61.9°), fcc-Co(002) (6'=25.76° and X'-54.8°) and 
ZnO(lOll) (6'=18.13° and X'^61.6°) reveal the in-plane 
orient relationship for Co NCs respect to ZnO. 
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tion response is detected for the sample annealed at 1073 
K. Note that increases drastically above 330 K af- 
ter annealing at 923 K. However Jacobsohn et al. [56] 
reported a much lower (250 K) of hep Co NCs with 
similar grain size (~10 nm). The higher could be 
due to some other anisotropy energy, which stablized the 
superparamagnetism at higher temperature. This will be 
discussed in section [Vl 

Hysteresis loops were measured for all samples in both 
parallel (ZnO[1010] ||hcp-Co [1010] ||fcc-Co [112]) and per- 
pendicular (ZnO[0001]||hcp-Co[0001]||fcc-Co[lll]) direc- 
tions. Figures flT[b)-(d) reveal that the anisotropy and 



FIG. 10: (a) 20-0 scans for ZnO(lOll) and hcp-Co(lOll), 
those small sharp peaks are from artificial noise in the mea- 
surements; (b) 20-0 scan for fcc-Co(002). 



coercivity can be tuned by different annealing proce- 
dures. The as-implanted sample only consists of hep Co 
NCs, which persist the bulk-like anisotropy along [0001] 
direction with a high coercivity (Figure [iTJb)). The sam- 
ple annealed at 823 K mainly consists of fcc-Co with easy 
axis is along the fcc-Co [112] direction (Figure [TTJc)). Af- 
ter annealing at 923 K, the sample shows an easy axis 
both along hcp-Co[1010] and hcp-Co[0001], but a higher 
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FIG. 11: (a) ZFC/FC magnetization curves at 50 Oe for the 
samples after implantation and annealing at different temper- 
atures. Solid symbols are FC curves, while open symbols are 
ZFC curves; (b)-(d) M-H curves measured at 5 K for all sam- 
ples: along ZnO[1010]||hcp-Co[1010]||fcc-Co[112] (solid sym- 
bols) and ZnO[0001]||hcp-Co[0001]||fcc-Co[lll] (open sym- 
bols). Implantation or annealing temperature is shown. 
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FIG. 12: (a) ZFC/FC magnetization curves at 50 Oe for the 
samples after implantation and annealing at different tem- 
peratures. Solid symbols are FC curves, while open sym- 
bols are ZFC curves; (b) M-H curves measured at 5 K for all 
samples along ZnO[1010] ||fcc-Ni[112]; (c) and (d) M-H curves 
measured at 5 K for all samples along ZnO[1010] ||fcc-Ni[112] 
(solid symbols) and ZnO[0001] ||fcc-Co[lll] (open symbols) af- 
ter 823 K and 923 K annealing, respectively. Implantation or 
annealing temperature is shown. 



coercivity along the latter direction. The magnetic prop- 
erties of different samples are listed in table [III 



Figure 12 show the magnetic properties of Ni im- 
planted ZnO upon annealing. As expected, the 823 K 
annealing increases in the ZFC curve due to the in- 
crease of crystallite size. However annealing at 923 K, 




5(26) (deg.) 

FIG. 13: XRD scan of Ni(200) and Ni(lll) 



the shape of ZFC magnetization curve is deviated from 
others. There is no real maximum, but a broad plateau. 
The formation of Ni/NiO core/shell structure could in- 
troduce the exchange coupling, which contributes an- 
other anisotropy energy. This will be discussed in the 
following section. Note that after annealing at 1073 K the 
sample still exhibits a non- neglectful response in ZFC/FC 
magnetization measurement. There are could be a small 
amount of Ni nanocrystals remaining, while they are be- 
yond the detection limit in SR-XRD measurement. 



V. DISCUSSION 
A. Magnetic anisotropy of Co and Ni NCs 

As shown before, the magnetic anisotropy of Ni NCs 
embedded in ZnO are drastically different from the bulk 
Ni. For bulk Ni, the magnetocrystalline anisotropy con- 
stant of Ki is -5.7x10^ J m-3 at 300 K. The [111] direc- 
tion is the easy axis. For a single magnetic NC, the mag- 
netocrystalline, shape, and magnetoelastic anisotropy 
have to be considered. 

In principle, a uniformly magnetized single domain 
spherical particle has no shape anisotropy, because the 
demagnetizing factors are isotropic in all directions. 
However, in the case of a nonspherical sample it will be 
easier to magnetize along a long axis than along a short 
one. The FWHM of the Ni peak in XRD scans is 
a measure of the crystallite size. Along the Ni[lll] di- 
rection, the crystallite size is estimated to be around 8 
nm. For the in-plane direction, the crystallite size could 
be estimated by measuring a diffraction plane not par- 
allel with sample surface. We chosen Ni(200) diffraction 
plane. Figure [13] shows the comparison of the normalized 
Ni(lll) and (200) peaks. Actually, the Ni(200) is broader 
than (111), which means the crystallite size along the in- 
plane direction is even smaller than the [111] direction. 
Therefore, the shape anisotropy is not the key reason to 
induce the easy axis along the in-plane direction. 

Now we consider the strain anisotropy. This kind of 
anisotropy is often described by a magnetoelastic energy 
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term 



B. Annealing at 923 K 



^magnetoelastic 2 COS 



(3) 



where is the magnetostriction constant, a the stress, 
and the angle between magnetization and the strain 
tensor axis [57^. The strain e along [111] direction of our 
Ni NCs can be calculated according to XRD measure- 
ment. 

The lattice spacing is given by Bragg law: 



2dexp sin t 



A 



(4) 



where dexp is the lattice spacing, the Bragg angle, and 
A=0.154 nm the wavelength of x-ray. is obtained from 
Figure [2] The strain is defined as following 



{dexp -d)/d 



(5) 



where d is the theoretical lattice spacing for bulk Ni or 
Co. 

Using the approach and parameters in Refs. [58*, W, 
the magnetoelastic anisotropy constant is calculated (see 
table III). For Ni NCs, the magnetoelastic anisotropy 



constant is one order higher than the magnetocrys- 
talline anisotropy constant. Therefore the magnetoelas- 
tic anisotropy energy dominates the total anisotropy en- 
ergy. This finding demonstrates the possibility to tune 
the magnetic properties by embedding magnetic NCs in 
different host matrix. The annealing at 823 K does not 
change the strain status significantly. After annealing up 
to 923 K, the elastic strain is partially released. However, 
the magnetoelastic anisotropy energy is still much higher 
than the magnetocrystalline one. Note the anisotropy 
change in Figure 12 ^d) after 923 K annealing by compar- 
ing with the 823 K annealing. There could be another 
anisotropy source after 923 K annealing, which will be 
discussed in section [VBl 

For hcp-Co, the magnetocrystalline anisotropy con- 
stant is very large, and it is difficult to be dominated by 
other anisotropy energy contribution (as shown in Table 
Therefore the hcp-Co NCs in the as-implanted state 
still persist the bulk like anisotropy behavior. However, 
the annealing at 823 K resulted in coexistence of fcc-Co 
and hcp-Co NCs. Since fcc-Co is a meta-stable state, 
its magnetic and mechanic data are very limited. There- 
fore in this paper, we will not discuss its magnetoelastic 
anisotropy energy. Nevertheless fcc-Co [112] is the easy 
axis rather than [111], which could be due the lattice 
strain, as the same reason for Ni nanocrystals in ZnO. 
After annealing at 923 K, hcp-Co is the only phase, and 
the elastic strain is not significantly changed compared 
with the as-implanted state (see Figure [sja)). However, 
the anisotropy is obviously changed. Like the Ni case, 
there could be another anisotropy source, which will be 
discussed in section IV Bl 



After annealing at 923 K, the magnetic proper- 
ties, namely ZFC/FC magnetization, and magnetic 
anisotropy, are changed significantly compared with 
other samples. For the Co case, Tb in the ZFC curve 
increases drastically, while the FC curve is not mono- 
decreased with increasing temperature. For Ni case, 
there is a broad plateau in ZFC curve. 

One explanation for the higher Tb is the increase and 
broadening of the Co crystallite size. Jacobson et al. cal- 
culated the ZFC curves by varying different parameters, 
including the size distribution [52 . It is found that a 
slight broadening can result in a very broad and high 
ZFC curve. This can well explain the ZFC magnetiza- 
tion, but not the FC magnetization. Moreover the dras- 
tic changes in coercivity, and magnetic anisotropy cannot 
be explained in such an approach. In addition. Table [ll| 
lists the comparison of the XRD peak area and the mag- 
netization upon annealing. Note that the XRD peak area 
decreases drastically by almost one order after 923 K an- 
nealing, which indicates the decreasing of the amount of 
the Co or Ni NCs, however the saturation magnetization 
only decreases slightly. 

Another explanation is due to other anisotropic en- 
ergy contribution. Skumryev et al. have found that the 
exchange coupling between Co NCs and their CoO shells 
drastically increase the blocking temperature [60 . In our 
XRD measurement, the amount of metallic Co and Ni de- 
creases after annealing at 923 K, which is very probably 
due to the oxidation. Therefore we could assume the for- 
mation of a Co/CoO core/shell structure in our sample, 
and the exchange coupling increases the blocking tem- 
perature. The exchange coupling is further confirmed by 
the vertical shift of the magnetization loop (Figure 14). 
The loop-curve is shifted along the magnetization axis 
after cooling from 350 K in an applied field. 

The vertical shift of the magnetization loop is a strong 
evidence for the presence of an interfacial interaction 
between an anitferromagnet (AFM) and a ferromagnet 
(FM) [6T1, i.e. the hep Co NCs are surrounded by CoO, 
which is amorphous, and therefore cannot be detected by 
XRD. Dobrynin et al. [61] presented a model to discuss 
the exchange coupling of nanoscaled Co/CoO core/shell 
structures. Below a critical size (12 nm) of Co cores, the 
interfacial exchange energy is larger than both the Zee- 
man energy of FM and the anisotropy energy of the AFM 
due to a large surface-to- volume ratio of NCs, and conse- 
quently some spins in the FM part can be frozen by the 
AFM part, leading to a vertical shift along the magneti- 
zation axis after field cooling. This model well explains 
our Co/ CoO system with an average diameter of 10 nm 
for Co NCs. The vertical shift of the magnetization loop 
decreases with increasing temperature and disappears at 
a temperature between 200 and 300 K. This is consistent 
with the Neel temperature of 290 K for CoO [60]. Such 
a vertical shift is also observed in the FC-loop measured 
along ZnO[1010]||hcp-Co[1010] (in-plane) (not shown). 
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TABLE III: Comparison of magnetic-anisotropy energy density: magnetocrystalline Ki, magnetoelastic Kme- As is the 
magnetostriction coefficient, E the Young modulus, e the elastic strain, and a the stress {a = Ee). 



Nano crystals 


kME (10^ Jm"^) 


ki (10^ Jm-^)^ 




E (GPa)" 


e 


a (Gpa) 


Ni (As-imp.) 


54 


-5.7 


-24x10"^ 


200 


0.011 


2.2 


Ni (823 K ann.) 


50 


-5.7 


-24x10"^ 


200 


0.010 


2.0 


Ni (923 K ann.) 


25 


-5.7 


-24x10"^ 


200 


0.005 


1.0 


hcp-Co (As-imp.) 


125 


500 


-5.5x10"^ 


209 


0.011 


2.3 



^Data from Refs. [58] [59] . 
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FIG. 14: The sample annealed at 923 K shows a vertical shift 
along the magnetization axis after field cooling (11=2000 Oe) 
due to the exchange coupling between FM and AFM materi- 
als. The field is along ZnO[0001] direction. The inset shows 
the temperature dependence of the vertical shift. 



For the Ni implanted sample annealed at 923 K, the 
hysteresis loop was also measured under both ZFC and 
FC conditions (not shown). A similar shift along the 
magnetization axis is observed. Therefore we would as- 
sume the formation of a Co/CoO (Ni/NiO) core/shell 
structure. The exchange coupling between Co (Ni) and 
CoO (NiO) contributes another anisotropy energy. This 
explains the high blocking temperature and the change 
in magnetic-anisotropy after annealing at 923 K. 

Note that for both of cases, the magnetization is only 
slightly decreased after annealing at 923 K. This is con- 
troversial with the drastic decrease of the diffraction in- 
tensity (see Figure [8| . One possible reason is the for- 
mation of small Co or Ni nanocrystals due to annealing. 
As shown in Table |lj only around 10-40% of Co or Ni 
form as metallic nanocrystals and others remain as dis- 
persed ones. Annealing at 923 K, on one hand, oxidized 
some metallic nanocrystals partially, on the other hand, 
could induce the gettering of dispersed Co and Ni and re- 
sult in small metallic nanocrystals. Below a critical size, 
nanocrystals are non-detectable by XRD. Another reason 
could be that only a thin shell of Co or Ni nanocrystals 
transforms into oxides, which results in the invisibility 



of oxides in XRD even they are crystalline. Obviously 
a detailed investigation using transmission electron mi- 
croscopy or other more sensitive techniques should be 
performed to clarify this controversy. 



VI. CONCLUSIONS 

A thorough characterization of the structural and mag- 
netic properties has been presented on Co and Ni im- 
planted ZnO single crystals. The results by SR-XRD and 
SQUID magnetometery correlate well with each other. 
The major conclusions are summarized as follows. 

(1) Co or Ni NCs have been formed in ZnO upon ion 
implantation. Their crystalline sizes, generally below 10 
nm, increase with increasing fluence. From several to 
37% percent of implanted Co or Ni is in metallic states, 
while the remaining could be diluted into ZnO matrix. 
The Co or Ni NCs are the origin of the measured ferro- 
magnetism. 

(2) The Co and Ni NCs are crystallograph- 
ically oriented with respect to ZnO host ma- 
trix. The orientation relationship is as fol- 
lows: hcp-Co(OOOl) [1010] II ZnO(OOOl) [1010] , and 
Ni(lll)[112]||ZnO(0001)[1010]. This well ordered 
structure of NCs could result in a rather smooth inter- 
face between them and ZnO host, and makes the hybrid 
of ferromagnetic NCs and semiconductors promising for 
spintronics functionality. 

(3) Magnetic anisotropy is observed for Co or Ni NCs in 
ZnO. Especially for the Ni NCs, the anisotropy is differ- 
ent from the bulk crystals. The extra anisotropy energy 
is attributed to the lattice strain impressed from the host 
matrix. This opens a route to artificially tune the mag- 
netic properties of nanoparticles by selecting of substrate 
materials. 

(4) The structure and magnetic properties of Co or 
Ni NCs embedded inside ZnO can be tuned by post- 
annealing. For the Co case, 823 K annealing results 
in the co-exists of fcc-Co and hcp-Co. The magnetic 
anisotropy is changed from out-of-plane to in-plane. An- 
nealing at 923 K could have partially oxidized metallic 
Co and Ni, and result in Co/CoO (Ni/NiO) core/shell 
structures. After annealing at 1073 K, no Co or Ni NCs 
can be detectable within the detection limit of SR-XRD, 
at the same time, the samples shows no pronounced fer- 
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romagnetism down to 5 K. 
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